Valence bond solid states with symplectic symmetry 
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We introduce a one-dimensional valence bond solid (VBS) state with symplectic symmetry SP(n) 
and construct the corresponding parent Hamiltonian. We argue that there is a gap in the spectrum. 
We calculate exactly the static correlation functions, which fall off exponentially. Hence the model 
introduced here shares all properties of the Haldane scenario for integer-spin quantum antiferromag- 
nets. We further show that the VBS state possesses string order and discuss its generalization to 
higher dimensions. 
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I. INTRODUCTION 

In 1987 Affleck, Kennedy, Lieb, and Tasaki (AKLT) [1] 
introduced the SU(2) valence bond solid (VBS) state and 
showed that it is the unique ground state of a special an- 
tiferromagnetic spin-1 chain. This model possesses all 
properties of the Haldane scenario for integer-spin quan- 
tum antiferromagnets 0], namely a unique ground state, 
an energy gap between the ground state and the exci- 
tations, and exponentially decaying correlations in the 
ground state. Soon after its discovery the VBS state was 
reformulated in terms of Schwinger bosons [1, 0] . This 
formulation revealed a striking analogy between the VBS 
state and the Laughlin state in the fractional quantum 
Hall effect Q and enabled the analysis of the excitations 
above the AKLT state using a single-mode approxima- 
tion Q. Following these developments the AKLT model 
was widely used to study general properties of spin-1 
chains, for example the appearance of hidden string or- 
der and a Z2 X Z2 -symmetry breaking in the Haldane 
phase This success has also motivated the study of 
g-deformed AKLT chains Q as well as SU(n) generaliza- 
tions of the VBS construction Very recently, Tu 
et at extended the investigation of hidden string order 
to SO(n) symmetric Hamiltonians [Tlj . 

Another important invention was the formulation of 
generalized VBS states in terms of finitely correlated or 
matrix product states [T^ . [isj , which in particular allows 
the relatively easy calculation of correlation functions. 
Ostlund and Rommer [13] showed that the wave func- 
tions appearing in the density matrix renormalization 
group (DMRG) method (Tsj are represented by matrix 
product states. Since then DMRG algorithms, which 
make direct use of the matrix product state formula- 
tion [l^, have been developed. 

On the other hand, large-n techniques based on 
symplectic symmetry were introduced by Read and 
Sachdev [l3] to study frustrated antiferromagnets on a 
square lattice. If one places symplectic spins transform- 
ing under a given representation of SP(n) on an arbi- 
trary lattice, it is always possible to form singlet bonds 



between any two sites. This is not true for unitary spins 
transforming under SU(n), where the formation of sin- 
glets is in general only possible on bonds between a rep- 
resentation and their complex conjugate representation. 
This restricts the applicability of SU(rt) techniques to 
bipartite lattices. The SP(n) technique was afterwards 
widely used to study frustrated antiferromagnets on var- 
ious lattices fisl , doped antiferromagnets [l^ , paired 
Fermi gases [20|, stripes in high-temperature supercon- 
ductors [2H, an d heavy- fermion systems [2^. Recently, 
Flint et al. |23| introduced the "symplectic-n" approach 
which links time reversal and symplectic symmetry of 
spins by eliminating unwanted dipole moment operators 
in the decoupling procedure. This enabled the treatment 
of superconductivity on an equal footing with the Kondo 
effect. Moreover, Wu et al. [2^] pointed out that the 
model of ultra-cold spin-3/2 fermions with contact inter- 
action enjoys a generic SP(4) symmetry, which lead to 
further applications of the symplectic symmetry in the 
context of ultra-cold fermionic gases (25j . 

In this paper we will combine these aspects and gener- 
alize the VBS state to symplectic symmetry. We derive 
an exact parent Hamiltonian and argue that there ex- 
ists a finite gap in the excitation spectrum. We then 
use the representation of the VBS state in terms of a 
matrix product state to calculate the static correlation 
functions and the expectation values of various string op- 
erators. Finally we discuss the VBS state and possible 
parent Hamiltonians on higher-dimensional lattices. 



II. SYMPLECTIC SYMMETRY 

One of the key features of the group SU(2) is that two 
spins of arbitrary length S can always combine into a 
singlet, which is an essential condition for a proper de- 
scription of frustrated antiferromagnetism. The analog 
statement is not true for spins transforming under SU(n) 
with n > 3, where one has to deal with the tensor prod- 
uct of a representation and its complex conjugate one in 
order to form a singlet. In the language of antiferromag- 



netism this requires a bipartite lattice structure where 
one can place spins transforming under one representa- 
tion of SU(n) on one sublattice and the complex conju- 
gated spins on the other sublattice. One way to overcome 
this problem JJ] is the generalization of SU(2) spins to 
spins transforming under the symplectic group SP(7i), for 
which the formation of a singlet from two spins is always 
possible. 

The symplectic group SP(n) is the set of all unitary 
n X n-matrices U such that [2a, [13] 

U^IU = I, (1) 

where * denotes the transposed matrix and 
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As the matrix / is built up from blocks of 2 x 2 matrices, 
n has to be even. The generators of SP(n), which we 
denote by A"", a = 1, . . . , n{n + l)/2, have to satisfy 

+ = 0. (3) 

The elements in the group are obtained by U — 
exp (i X^a^a^") with real parameters 9a- The matrices 
A" play the same role as the Pauli matrices for SU(2) 
and equal them in the case n — 2. Hence there exists 
an isomorphism between SP(2) and SU(2); in particular 
the representations of SP(2) equal those of SU(2). An ex- 
plicit representation of the matrices A"" for SP(4) is given 
in App. [XI The irreducible representations of SP(n) can 
be labeled [1^ by (Ai . . . Xn/2), where the non-negative 
integers Ai have to satisfy Ai > A2 > • • ■ > A„/2- Explicit 
formulas for the dimensions of the irreducible representa- 
tions, the eigenvalues of the quadratic Casimir operator, 
and results on the decomposition of tensor products of ir- 
reducible representations are stated in App. [B) In Tab. |T] 
we have tabulated these properties for those irreducible 
representations which we will use to construct the VBS 
chain below. 

As a side note we mention that the symplectic g roup 
SP(n) naturally arises in Hamiltonian mechanics [28j . 
The n-dimensional phase space M contains the gener- 
alized coordinates qi, . . . , qn/2 ^nd their conjugated mo- 
menta pi, . . . ,Pn/2, which implies that n has to be even. 
The Hamiltonian : M — > R induces the time evolution 
via its vector field. The phase space is equipped with a 
skew-scalar product on its cotangent bundle T*M, i.e., 
a bilinear map (., .) : T*M x T*M R which satis- 
fies {x,y) = —{y,x). This skew-scalar product defines 
a volume element on the phase space. The symplectic 
group is now the set of all linear transformation under 
which this skew-scalar product is invariant. In particu- 
lar, the time evolution generated by the Hamiltonian is 
a symplectic transformation, which implies for example 
Liouville's theorem. 
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TABLE I; Simplest irreducible representations of SP(n), their 
dimensions, and the eigenvalues of the quadratic Casimir 
operator J^. We note that the representations (110... 0), 
(220 ... 0), and (310 ... 0) do not exist for SP(2) ^ SU(2). In 
this case the remaining representations are the singlet 0, the 
spinor representation i, the triplet 1, and the spin-2 repre- 
sentation 2. 
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FIG. 1: Weight diagram of the fundamental representation 
of SP(4). The states are labeled using the bosonic creation 
operators introduced in Q. and denote the diagonal 
spin operators, their eigenvalues are easily obtained using ((6|. 

III. DIMER CHAIN 

As a warmup exercise we first construct the SP(n) gen- 
eralization of the Majumdar- Ghosh model Let us 
consider a chain with N lattice sites and periodic bound- 
ary conditions, where we assume N to be even. On each 
lattice site we place an SP(n) spin transforming under 
the fundamental, n-dimensional representation (10 ... 0). 
A basis at each lattice site i may be written in terms of 
bosonic creation and annihilation operators ^ and foo-,z 
as [13 

k>. = <»|0>,, fT = l,...,n, (4) 

where |0)j denotes the vacuum at site i. The weight 
diagram of the fundamental representation of SP(4) is 
shown in Fig.[T] The action of the SP(n) spin operators 
J" on these basis states is given by 

1 " 

J^=2 E bUA^,\,,, (5) 

(7,(T' — 1 
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where we have introduced the vector notation J = 
(J\..., J"("+i)/2)'. The eigenvalue of the quadratic 
Casimir operator on each lattice site equals jf — {n + 
l)/4. Using the explicit expressions for the generators of 
SP(4) given in App. [Xlone finds for example 
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(6) 

We stress that the operators J^, J^, and as well as 
J®, J®, and J^" span two copies of the usual spin algebra 
su(2). 

On this SP(n) chain we consider the two linearly inde- 
pendent states represented by 

I o — o o — o o — o o — o ) "odd" 

(7) 

I o o — o o — o o — o o ) "even" 

where the symbol o — o stands for an SP(n) singlet or 
dimer formed by the spins on two neighboring lattice 
sites. In the state labeled as "odd" the SP(n) singlets 
are located on the bonds (2z — 1, 2i), whereas in the state 
labeled as "even" the SP(n) singlets are located on the 
bonds {2i,2i + 1). In the "even" state the right- and 
left-most spins also form an SP(n) singlet due to the as- 
sumed periodic boundary conditions. In order to con- 
struct a parent Hamiltonian, i.e., a Hamiltonian which 
has the two states © as its unique ground states, we 
note that the total SP(n) spin on each three neighboring 
sites has to contain a singlet and thus transforms under 
the fundamental representation (10. . .0). Hence, for all 
lattice sites i the operator {Ji + Ji+i + Ji+2)'^ — (n+l)/4 
annihilates the dimer states ([7]), and by taking the sum 
over all lattice sites we arrive at 

ffdimcr = iJ i+1 + -^JiJi+2 H ^ ^ ■ (8) 

We have checked numerically for n = 4 and = 8 that 
the model ([5]) possesses exactly two zero-energy ground 
states. For n = 2 one obtains the original Majumdar- 
Ghosh model 



IV. VBS CHAIN 

In this section we construct the SP(7i) VBS state on a 
chain and derive the corresponding parent Hamiltonian. 
In the next sections we will then discuss the excitations 
above the VBS state, its static correlation functions, and 
the appearance of string order. 

Let us consider again a chain with A^ lattice sites and 
periodic boundary conditions, but now A^ may be even 
or odd. At each lattice site we place two copies of the 
fundamental representation (10 ... 0), i.e., we obtain the 
tensor product (the decomposition of tensor products in 
irreducible representations was derived in Refs. (sil . [3^ 
and is presented in App. [B| 

(10 ... 0) (10 ... 0) = (20 ... 0) e (no ... 0) e (0 ... 0). 

(9) 

We note that for n = 2 the representation (110 ... 0) does 
not exist and we recover ^(g)i = l©0. In the tensor 
product © we now project onto the adjoint, n{n + l)/2- 
dimensional representation (20... 0). An explicit basis 
for this representation can be constructed [30] from the 
bosonic basis of the fundamental representation For 
n — 4 this basis will be stated explicitly in Sec. lVIl With 
this procedure we construct a chain of adjoint representa- 
tions, which is the direct generalization of a spin-1 chain 
for SU(2). If we consider the total SP(n) spin of two 
neighboring sites we find the decomposition 

(20 ... 0) ® (20 ... 0) = (40 ... 0) e (310 ... 0) 

® (220 ... 0) ® (20 ... 0) © (110 ... 0) ® (0 ... 0). 

(10) 

For n — 2 the second, third, and fifth representation on 
the right-hand side do not exist and pO)) simplifies to 
1(8)1 = 2010 0. 

Starting with such a chain of adjoint representations, 
we can construct the VBS state as follows: We form a 
singlet between one of the fundamental representations 
(10 ... 0) on lattice site i with one of the (10 . . . 0)'s on the 
neighboring site i — 1 while we form another singlet with 
the second representation (10 ... 0) on lattice site i with 
one of the (10 . . .0)'s on the neighboring site i -|- 1. We 
stress that the formation of these singlets is imposed in 
addition to the already implemented projection onto the 
adjoint representation at each lattice site. If we further 
impose periodic boundary conditions this yields a unique 
VBS state |^'vbs)j which is translationally invariant and 
can be represented graphically as shown in Fig. O 

The parent Hamiltonian for the VBS state is con- 
structed by noting that on each two neighboring sites 
in the VBS state we find one singlet and two uncoupled 
fundamental representations. Hence, the total SP(7i) 
spin on two neighboring sites is given by the tensor 
product (10 . . . 0) «) (10 . . . 0) given in If we con- 

struct an operator which is identical to zero on ^ but 
takes strictly positive values on the complement of ^ 
in (20 ... 0) (g) (20 ... 0), we will obtain the VBS state as 
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FIG. 2: Graphical representation of the VBS state |^vbs), 
the unique ground state of (|12|) . Each circle stands for a 
fundamental representation (10...0), each line joining two 
circles for a singlet bond, and each oval for a lattice site on 
which we project onto the adjoint representation (20 ... 0). 



zero-energy ground state. This operation is most eas- 
ily implemented using the quadratic Casimir operator 
(Ji + Ji+i)^ on the bond + which takes the values 
(n -I- 2)/2, n/2, and on the representations in ^ and 
n + 4, n + 2, and n + 1 on the remaining representations 
in (|10p , respectively. Explicitly we will use on each bond 
(i,i + l): 



Pi,i+1 — 



1 



n 5n2 -I- 26n + 32 



J, + J,;+l 



J, 



(11) 



We stress that the operators are not simple pro- 

jectors, as Pi,i+i takes different values on the subspaces 
(220... 0), (310... 0), and (40... 0). We note that for 
n = 2 the last factor in (fTT|) is not necessary as the cor- 
responding representation (110 ... 0) does not exist. The 
normalization of Pi^i+i is chosen in order to obtain a fi- 
nite expectation value for the energy of each individual 
bond in the limit n — > oo. In this limit the operator pi|) 
becomes an orthogonal projector (up to the multiplica- 
tive constant 1/10) onto the complement of ([9]) in (fTO| . 
The parent Hamiltonian for the VBS state I^E'vbs) is now 
obtained hy H — Pi^i+i together with = (n + 2)/2: 
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Here the operators Jf live in the adjoined representa- 
tion and can be represented by n(n -I- 1)/2 x n{n + l)/2- 
matrices. As the operator ([TT|) takes strictly positive 
values on (220 ... 0), (310 ... 0), and (40 ... 0), all states 
except the VBS state are lifted to higher energies. We 
have checked numerically for n — A and = 3 that the 
VBS state is the unique ground state of p^ . A proof of 
the uniqueness can be obtained by generalizing the proof 
of the uniqueness of the ground state of the q-deformed 
VBS model Q. The Hamiltonian contains cubic terms 
as we had to use three factors in the operators (|lip . As 
explained above the third factor is superfluous for n = 2, 
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FIG. 3: If two spinous (represented by the full circles) move 
apart from each other, the region between them will not have 
the same structure as in the ground state (we have sketched 
one way of how the SP(n) spins may rearrange into singlets). 
This causes an energy cost proportional to their distance and 
results in spinon confinement. 



omitting it yields the orig inal AKLT model fi\. By keep- 
ing the third factor, however, we obtain an alternative 
parent Hamiltonian for the spin-1 VBS state. 

The VBS construction described above can also be 
done for a chain with open boundary conditions. In this 
case we are left with one uncoupled fundamental repre- 
sentation at each end of the chain and we hence find r? 
linearly independent VBS states. The parent Hamilto- 
nian for these states is given by with the summation 
restricted tol<i<A^— 1. 



EXCITATIONS AND ENERGY GAP 



The Hamiltonian p2)l was constructed to be the exact 
parent Hamiltonian for the VBS state I^'vbs)- Although 
its ground state is known in all detail it is much harder to 
get results on the excitations above it. The simplest oper- 
ation on the state I^'vbs) one can imagine is to break one 
of the singlets, say the singlet on the bond {i, i+1). Doing 
so we find two uncoupled SP(n) spins each transforming 
under the fundamental representation, which we will call 
spinous in the following. The resulting state is clearly not 
an eigenstate of (HH). Nevertheless, the spinous are use- 
ful to perform the following Gedankenexperiment: Let us 
pin the first spinon at site i and move the other spinon 
to the right (see Fig. [3]). The region between them has 
now a different structure than the ground state and is 
not annihilated by As the energy cost grows lin- 

early with the distance, the spinons are subject to a lin- 
ear confinement potential and hence can only appear in 
bound states. The relative motion of the spinons will be 
described by a non-harmonic oscillator whose zero-point 
energy yields a finite gap for the creation of spinon-spinon 
bound states. This is consistent with the picture that the 
origin of the Haldane gap is a confinement force between 
spinons E^, [s^ . A similar argumentation was applied by 
Greiter [3^ to the excitations of the two-leg t-J ladder. 
Although this Gedankenexperiment suggests the appear- 
ance of an energy gap, we stress that the spinon bound 
states may not constitute good trial wave functions for 
the actual low-lying excitations in the model. 

A possible way to prove the existence of a gap above 
the ground state is provided by the extension of results 
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by Knabe [3^ on a class of SU(2) VBS Hamiltonians in- 
cluding the original AKLT model (details of the deriva- 
tion are given in App.lCjl. Let us consider a Hamiltonian 
of the form 



H 



N 

E 



(13) 



with periodic boundary conditions. We assume that < 
Pi.i+i < 1 as well as the existence of at least one zero- 
energy ground state of (llSp . The idea is to establish the 
inequality 



>eH, e> 0, 



(14) 



which implies that the lowest non- vanishing eigenvalue of 
H is larger than e. As we show in App. [Cl Eq. (|14p can 
be derived if the same model on a chain with m + 1 sites 
and open boundary conditions satisfies 



> 



1 



where m > 2 and 



z+m — 1 
k—i 



k,k-\-l- 



(15) 



(16) 



Hence the proof is finished if we can show that is 
satisfied for a suitable integer m. This was achieved in 
Ref. m for SU(2) VBS chains with spins S = 1,3/2,2, 
and 5/2 on each lattice site. Unfortunately, exact diag- 
onalization of the SP(4) model with open boundary 
conditions for m = 2 suggests that in order to establish 
the inequality > 1 /m one has to study chains with at 
least ten lattice sites. 



VI. STATIC CORRELATION FUNCTIONS 

The VBS state 1^1/ vbs) can be written as a matrix prod- 
uct state. We will restrict ourselves to the case n = 4 in 
the following. A suitable basis for the adjoint represen- 
tation of SP(4) at lattice site i can be obtained from 
the bosonic basis of the fundamental representation in- 
troduced above as [301 
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FIG. 4: Weight diagram of the adjoint representation of 
SP(4). The state with = J^" = is doubly degenerate. 
The states are labeled using the bosonic creation operators as 
in ini). 



In order to derive the matrix product representation 
we first rewrite the singlet on the bond + 1) as 



"l.i"2.i+l "2.i"l. 



"3,i"4,i+l ".1 



4,1^^3, i+l 



b'^ b'^ b'^ b^ 



(18) 



Second, at each lattice site i we use the outer product to 
combine the two vectors originating from the rewriting 
(|18p on the bonds (i — l,i) and + 1) into a matrix 
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(19) 

Assuming periodic boundary conditions the VBS state 
can then be written as the trace of the matrix product 



N 



|*VBs) = trmM, 



(20) 



We have illustrated these basis states in the weight dia- 
gram of the adjoint representation shown in Fig. 3) 



Starting from this representation the static correlation 
functions in the SP(4) VBS state can be calculated by 
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applying the method introduced by Kliimper et al. for 
the analysis of the g-deformed model [8i] . As the first step 
we calculate the norm of the VBS state. This is done by 
introducing the complex conjugated matrix M accord- 
ing to Ma-a' = M*^,, i.e., by simply taking the complex 
conjugate of each matrix element in (fT9|) without trans- 
posing the matrix. We then define the f 6 x 16 transfer 
matrix R at any lattice site as 



{aT).(a'T') 



M , M , 



(21) 



order was further associated with the breaking of a Z2 x 
Z2-symmetry in the Haldane phase and the appearance of 
a four- fold degenerate ground state on the open chain [3] ■ 
We will find a similar behavior in the SP(n) VBS model. 

In analogy to Refs. @ we define the string operators 



-Jf exp i^^5]j,M J]- 



(28) 



k=2 c 



where we order the indices as a,/3 = 1,...,16 ^ 
(11), (12), . . . , (44). The norm of the VBS state is now 
given by 



(^VBsl *VBs) = ^^[R'') = 5'' + 10 (-1)^^ + 5, (22) 

where we have evaluated the trace by diagonalization of 
R. In the second step we calculate the expectation value 
(*VBs| ^1 Jj l^'vBs)- We introduce the transfer-matrix 
representation of the spin operators by 

Jal3 — J(aT),(<T'T') = Mca' M^t' ■ (23) 

Here the operator acts on the elements of M as 

(24) 

which implies for example J^b\b2 = 0. This yields 



^-VBsl Jf J- I^VBs) = tr Ji R1-' J, R 



(25) 



which is easily evaluated by diagonalization of R. As the 
state l^'vBs) enjoys full SP(4) symmetry we arrive at 



6\ _ (*VBS| J^I^VBS) 
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_9_ + 
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(26) 



1+ 5^(2(-l)^^ + l) 



In the general case of SP(n) the same steps yield in the 
thermodynamic limit N 00 



ja jl 



5ab 



(n + l) 



(27) 



Here the correlation length is given by f = l/ln(n + 1) 
and vanishes in the limit n ~^ 00. We also recover the 
known result for the AKLT chain [l| . 



VII. STRING ORDER 

It is well known |6| that there exists a hidden nonlocal 
topological order or string order in the AKLT model. In 
fact, this order is found in the whole Haldane phase in 
the phase diagram of the general spin-1 chain. This string 



where the second sum is over all c for which is diagonal 
and Jf and Jj have to be diagonal as well. In the SP(4) 
model the summation is over c = 3, 10 and we have a, 6 £ 
{3, 10}. Using the transfer-matrix technique we obtain in 
the thermodynamic limit N 00 
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(29) 



which remain finite for arbitrary large values of j. In 
particular, the sum over all four expectation values ((29|) 
is independent of j. In analogy to the original AKLT 
model we expect this hidden string order as well as the 
16-fold degeneracy of the ground state of a chain with 
open boundary conditions to be a consequence of the 
breaking of a discrete symmetry (Z4 x Z4). 

We have also calculated the expectation values of the 
nine string operators (|28ll in the SP(6) model. Together 
with ^ and the result [gj] for SU(2) this leads us to the 
conjecture for general n: 
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, a^b. 



(30) 



Although each of the expectation values ([50]) vanishes 
in the limit n 00, the number of string operators in- 
creases and one obtains 
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a^b 



^ab 



1 

4' 



(31) 



where the sum is over all a and h for which J" and 
are diagonal. We note that f3T|) can be written elegantly 
as a single string operator by replacing Jf and jj* in ((28|) 
by the sum over all diagonal generators ^1 respec- 
tively. 
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singlet 
bond 



projection 
onto (30 . . . 0) 



FIG. 5: Graphical representation of the VBS state on a hexag- 
onal lattice. Each small circle represents a fundamental rep- 
resentation (10 ... 0), each line joining two circles for a singlet 
bond, and each large circle a lattice site on which we project 
onto the representation (30 ... 0). 



VIII. TWO-DIMENSIONAL VBS MODEL 

Finally we would like to discuss the VBS construction 
on higher-dimensional lattices. The simplest example is 
provided by the honeycomb lattice (with coordination 
number z ~ 3) with representations (30 ... 0) on each 
lattice site. The VBS state is obtained by placing three 
fundamental representations on each lattice site and pro- 
jecting onto the representation (30 ... 0) as well as form- 
ing a singlet of each one of them with a fundamental 
representation on a neighboring site (see Fig. [S]). Hence, 
on each bond we obtain a singlet formed in this way and 
four uncoupled fundamental representations. The corre- 
sponding tensor product decomposes as 

(10... 0)®^ = (40... 0)® 3 - (310... 0) 

® 2 • (220 ... 0) ® 3 • (2110 ... 0) ® 6 • (20 ... 0) 
® (11110...0)®6- (110...0)®3 - (0...0). 

(32) 

Since this tensor product contains in general eight dif- 
ferent irreducible representations, the "projection" oper- 
ator analog to pT|) and hence the Hamiltonian contains 

the bond operators (Ji Jj)'' with a power up to fc = 8. 
For the cases n = 2,4, and 6, however, some represen- 
tations on the right-hand side of (I32|) do not exists and 
one obtains powers up to A: = 3, 6, and 7, respectively. 
The explicit construction of the analog SU(2) model with 
5 = 3/2 spins on the lattice sites can be found in Refs. 
It is clear from the arguments above that the VBS con- 
struction on lattices with larger coordination number or 
in higher dimensions will lead to a parent Hamiltonian 
which contains even higher powers of the bond operators 



IX. CONCLUSIONS 

In conclusion, we have introduced a spin chain with 
symplectic symmetry SP(n) which shares all properties 
of the Haldane scenario for integer-spin quantum anti- 
ferromagnets: (i) a unique ground state, (ii) a finite gap 
in the energy spectrum above the ground state, and (iii) 
ground-state correlation functions which fall off exponen- 
tially. Furthermore we have shown that the ground state 
possesses string order. We point out that in the limit 
n oo the string order remains finite and the correla- 
tion length vanishes. The application of the large-n ap- 
proach to the considered models might be an interesting 
extension of this work. 
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APPENDIX A: EXPLICIT REPRESENTATIONS 
FOR SP(4) 

An explicit representation of the generators of SP(4) 
is provided by 



A' = 
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The normalization is chosen to be 

tr {A- A') = 2Sab 



(Al) 



The matrices v4°, a ~ 1, . . . , 10, form a basis of sp(4), 
the Lie algebra of SP(4). They satisfy the commutation 
relations 



[A'',A''] = 2/"''"^'= 



(A2) 



The structure constants f^^^ are totally antisymmetric 
and obey Jacobi's identity 



j^abc jcde jbac j 



bdc fcae 



dac fche 



0. 



(A3) 



Explicitly, all 84 non-vanishing structure constants are 
obtained by permutations of the indices from 



^123^^89,10^. 

yl56 _ j345 ^ j45,10 ^ j478 ^ j568 _ j579 _ j67,10 
jl47 _ y246 _ y257 _ j367 _ j469 _ _^ 



1 

2' 



(A4) 



sp(4) has rank two, the Cartan subalgebra is spanned by 
A^ and A^° . We note that sp(2)^ su(2) and sp(4)^ so(5). 
A possible matrix representation of the spin operators 
in the adjoint representation is given by (■Z");,^ = 
However, we stress that these matrices are not the rep- 
resentation matrices in the bosonic basis (|17p . 



APPENDIX B: SOME RESULTS ON THE 
REPRESENTATION THEORY OF SP(n) 

In this appendix we review some results on the rep- 
resentation theory of SP(n). First, the dimension of the 
irreducible representation (Ai . . . Xn/2) is given by the for- 
mula 

n/2 



dim[(Ai...A„/2)]=n'^^^^ 



=1 2 



i + l 



n/2 (Bl) 

nXi~\j-\-j — i\i+\j+n + 2 — i — j 



i<3 



n H- 2 — i — j 



Second, the eigenvalues of the quadratic Casimir oper- 
ator were derived by Nwachuku and Rashid and 
read using our conventions 



i=~n/2 * 2 2 j = -,i/2 ^ 



(B2) 



where for 1 < i < n/2 



+ i + A„/2+i- 



n — Ki. 



(B3) 



The following special cases allow a closed expression: 



{M...0): -{n + iy), 
uk , 



(B4) 



(;y_^0...0): — + + (B5) 

h terms 

The relation to Refs. [3^ is obtained by replacing n 
n/2 and rescaling the Casimir operator by a factor of 1 /8. 
For the simplest irreducible representations the formulas 
(lBT|) ~ ((B5t easily yield the results stated in Tab. [J 

Finally, we make use of the following results on the 
decomposition of tensor products into irreducible repre- 
sentations, which is in its general form due to Littel- 
mann |31 | and was specialized to the case we use here by 
Leung [32| : 



{pl^J.2 ■ ■ • Ain/2) (E> {v0...o) 



(B6) 



where the sum is over all integers ki, . . . , k„ subject to 
the conditions: 

Ki + . . . + Kn ^ I', 

< Ki < lli-i ~ — Kn^i+2 + Kn-i+1, 
< Kn-j < ^J'j + l — 
< K„/2+l < f-ln/2- 

where « = 2, 3, . . . , n/2 and j = 0, 1, . . . , n/2 — 2. 

APPENDIX C: DERIVATION OF EQ. (HH) 

In this appendix we will generalize results obtained by 
Knabe ^] on the existence of a gap in SU(2) VBS chains 
with arbitrary spin. Similar techniques were also used by 
Fannes et al. tl2| . The main difference of our result as 
compared to Ref. [35| is that the operators Pi^i+i are not 



assumed to be simple projectors. 

Let us start with (fT51) . The assumption Pi^i+i < 1 
yields -P^pi+i < -Pi,i+i, where inequalities between opera- 
tors are understood in the sense 



i+l 



< (V'IPm+iIV') < (V-l^) (Ci) 



for all states |^). In fact, the most useful results will 
be obtained if the largest eigenvalue of Pi^i+i equals one, 
which is obtained by multiplication of (llip with a suit- 
able constant. Using the definitions and one 
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easily finds 

= ^~ T ^li+l + ~; 7 ^i,m 

m — 1 ^ — ' TO — 1 -"^ — ' 

i=l i=l 
N 1 "'^ i+m-1 

+ X! + 1 ^ ^„ ;l X! X! Pk,k+lPl,l+l 

— l i—1 k,l—i 

\i-j\>l |fc-/|>l 

(C2) 

We can now use P^i+i < Pi,i+i, which impUes 



Y^i ^ ^1^^^ together with the fact that 

each of the terms Pi i^iPj j^i appears more often in the 
third sum than in the fourth sum. Therefore we get the 
inequahty 

H'>^y^h^^-^H. (C3) 

TO - 1 ^ TO - 1 ^ ' 

i=l 

Finally we can use as well as '^^hi^^ = mH to 
obtain m with e = ^ (e,„ ~ i) . 
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